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• • ■ ABSTRACT 
> ■ 

' Context. The earliest phases of massive star formation are currently much debated. 

H ' 

' Aims. In an effort to make progress, we took a census of Class 0-like protostellar dense 

cores in the NGC 3576 region, one of the nearest and most luminous embedded sites 
of high-mass star formation in the Galaxy. 

Methods. We used the P-ArTeMiS bolometer camera on the APEX telescope to produce 
the first 450 fim dust continuum map of the filamentary dense clump associated with 
NGC 3576. 

Results. Combining our 450 fim observations with existing data at other wavelengths, 
we have identified seven massive protostellar sources along the NGC 3576 filament and 
placed them in the Menv-Lboi evolutionary diagram for protostars. 
Conclusions. Comparison with theoretical evolutionary tracks suggests that these seven 
protostellar sources will evolve into massive stars with masses A/* ~ 15-50 Mq . Four 



* This publication is based on data acquired with the Atacama Pathfinder Experiment (APEX) 
in ESO program 080.C-0722. APEX is a collaboration between the Max-Planck-Institut fiir 
Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory. 
** Tables [T]& [5] and Figs. [2]&|3]are only available in electronic form. 
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sources are classified as candidate high-mass Class objects, two sources as massive 
Class I objects, and one source appears to be at an intermediate stage. 

Key words, stars: formation - stars: circumstellar matter - ISM: clouds - ISM: structure 
- ISM: individual objects (NGC 3576) - submillimeter 

1. Introduction 

The process of massive star formation remains poorly understood and is the subject of ma- 
jor ongoing debate (Zinnecker & Yorke 2007). Prom the observational perspective, there 
is a growing body of evidence to suggest that most, if not all, high-mass (M* > 8Mq) 
stars form in a clustered environment ( de Wit et al. 2005). There is no widely accepted 
evolutionary scenario for high-mass star formation, although a tentative sequence of empir- 
ical phases has been identified: from infrared dark cloud cores, to hot cores, hypercompact 
to ultracompact, compact and classical HIT regions. On the theoretical front, it is unclear 
whether high-mass stars form by direct cloud core collapse and subsequent accretion as a 
result of a scaled-up version of the process believed to be at work in low-mass protostars 
(McKee & Tan 2003), or by other processes such as competitive accretion and coalescence 
of intermediate-mass protostars in dense protoclusters (Bonnell et al. 2004). 

Decisive progress in this field will require high-resolution kinematical and modeling 
studies of high-mass protostellar objects and their embedded, cluster-forming environments 
at (sub)millimeter wavelengths. While such comprehensive studies are becoming feasible in 
the intermediate-mass regime (cf. Peretto et al. 2006, 2007), they are virtually impossible 
in the high-mass regime due to the lack of angular resolution and mapping speed of current 
submillimeter facilities. A major step forward will be provided by the upcoming ALMA in- 
terferometer when it becomes fully operational around 2013. In order to prepare for future 
ALMA studies, it is important to identify a number of candidate massive protostars and 
embedded protoclusters in the southern hemisphere. While several SIMBA 1.2 mm contin- 
uum surveys are available (e.g. Hill et al. 2005, Beltran et al. 2006), current searches for 
massive protostars in the southern sky lack both the angular resolution and submillimeter 
spectral coverage to provide an adequate characterization of the detected sources. 

At a distance of 2.4 kpc (Persi et al. 1994), NGC 3576 (also known as RCW 57 and 
G291.27-0.70) is one of the most luminous giant HII regions in our Galaxy at optical, 
infrared, and radio wavelengths (e.g. De Pree et al. 1999). It is embedded at the center 
of a filamentary dense clump observed in both 1.2 mm dust continuum emission (Hill et 
al. 2005) and tracers of dense molecular gas such as NH3 (Purcell et al. 2008). A deeply 
embedded near-IR cluster of more than 130 young stellar objects (YSOs) is known to be 
associated with this region (Persi et al. 1994). In addition, two 6.7 GHz methanol masers 
and a total of nine 22 GHz water masers (Caswell 2004, Purcell et al. 2008) have been 
detected along the extent of the cluster-forming clump. The NGC 3576 region is therefore 
an excellent laboratory for studying the embedded phases of massive star formation. Due to 
its southern declination, NGC 3576 is a potential target of choice for future high-resolution 
investigations with ALMA. 

Send offprint requests to: pandre@cea.fr 
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Here, we report the results of the first 450 /im dust continuum mapping observations 
of the NGC 3576 star-forming region, taken with the P-ArTcMiS bolometer camera on 
the APEX 12-m telescope. Section [5] describes the instrument and provides details about 
the observing run and data reduction. Section [3] presents our mapping results which are 
discussed in Section H) 

2. Observations and data reduction 

Our 450 /im observations of NGC 3576 were made in March 2007 and November 2007 with 
the P-ArTcMiS bolometer array on the Atacama Pathfinder Experiment (APEX) telescope 
located at an altitude of 5100 m at Llano de Chajnantor in Chile. P-ArTeMiS is a 16 x 16- 
pixel prototype of the ArTeMiS bolometer camera currently being built by CEA Saclay for 
APEX. ArTcMi^ is a large-format camera with a total of 5760 bolometric pixels which 
will operate simultaneously at 450 /im, 350 /im, and 200 /im in the Cassegrain cabin of 
APEX (e.g. Talvard et al. 2006, 2008). It uses the technology successfully developed by 
CEA for the PACS photometer instrument in the 60-210 fj,m wavelength regime on the 
Herschel Space Observatory (e.g. Simoens et al. 2004, Billot et al. 2006). Unlike classical 
bolometer arrays, such as LABOCA on APEX, the PACS and (P-)ArTeMiS instruments 
do not use feedhorns to concentrate the incoming submillimeter radiation, but act more 
like CCD cameras in the optical domain. The advantage of such a fiUed-array design is that 
the entire field of view of the camera is instantaneously fully sampled, and thus "jiggling" 
is not required to produce fully-sampled images. The P- ArTeMiS prototype was built with 
the objective of demonstrating the feasibility and potential of the full ArTeMiS instrument. 
P-ArTeMiS was first tested at 450 /im on the KOSMA 3m telescope in March 2006 and 
January 2007. It was then taken to Chile for initial tests on APEX during the second half 
of March 2007 under technical time allocated by the APEX Board, and then during ESO 
and Onsala scientific time in November 2007. P- ArTeMiS uses a detector array produced 
as part of the Herschel-PACS project, but not selected for the space mission, which was 
modified to operate at 450 /xm by adding a silicon-based filtering system on its surface 
(cf. Reveret et al. 2006). The P- ArTeMiS bolometers are operated at 300 mK inside a 
conventional He-refilling cryostat. During the 2007 runs, P- ArTeMiS was installed by our 
team at the visitor focus in the Nasmyth "A" cabin of the APEX telescope. 

Eight individual maps, corresponding to a total effective integration time of 3.6 hr, were 
obtained with P- ArTeMiS at 450 fim toward the NGC 3576 region using a total-power, on- 
the-fiy scanning mode. Each of these maps consisted of a series of scans in either Azimuth 
or Right Ascension taken with a scanning speed of 16"/sec or 20"/sec. The cross-scan 
step between consecutive scans was either 2" or 4". The sizes of the maps ranged from 
4' X 4' to 6' X 5'. The atmospheric opacity at zenith was monitored by taking sky dips with 
P- ArTeMiS and was found to be between 0.6 and 0.8 at A = 450 /im, corresponding to an 
amount of precipitable water vapor (PWV) between ^ 0.5 mm and ^ 0.7 mm. No dedicated 

^ See |http: / /www.apex-telescope.org/instruments/pi/artemis/ 1 
ArTeMiS stands for "ARchitectures de bolometres pour des TElescopes a grand champ de vue 
dans le domaine sub-MIUimetrique au Sol" in French. 
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pointing model was derived for P-ArTeMiS, but the pointing corrections proved to be very 
similar to, i.e. within ^ 10" of, the corrections appropriate for the APEX-2A instrument. 
Pointing, focus, and calibration measurements were achieved by taking both short 'spiral' 
scans and longer on-the-fly beam maps of Mars and Saturn. We estimated the absolute 
pointing accuracy to be ~ 5" and the absolute calibration uncertainty to be ^ 30%. 
The noise equivalent flux density (NEFD) per detector was 2 Jy.s^/^ . and the pixel 
separation between detectors on the sky was ^ 6". The individual maps of NGC 3576 were 
registered with a relative positional accuracy of ^ 2" and a relative calibration uncertainty 
of ^ 10%. As estimated from our maps of Mars (whose angular diameter was ~ 14" in 
November 2007), the main beam had a full width at half maximum (FWHM) ~ 10" and 
contained ^ 60% of the power, the rest being distributed in an "error beam" extending up 
to an angular radius of ^ 80". On-line data reduction was performed with the BoA software 
developed for LABOCA (SchuUer, Beelen et al., in prep.). Off-line data reduction, including 
baseline subtraction, removal of correlated skynoise and 1// noise, as well as subtraction 
of uncorrelated 1// noise using a method exploiting the high level of redundancy in the 
data (Roussel et al., in prep.), was performed with in-housc IDL routines. 



3. Mapping results and analysis 

Co-adding the eight individual P-ArTeMiS maps of NGC 3576, we obtained the 450 /im 
mosaic shown in Fig. [TJ At least seven localized dust emission peaks can be identified in 
this mosaic, which have been labeled here following the nomenclature of Purcell et al. 
(2008). They are distributed in four ^ parsec-scale clumps labeled SI, S3, S4, and S5. The 
positions and flux densities measured for these compact P-ArTeMiS sources are given in 
onhne Table H 

While the 450 /im emission peaks are well correlated with the 1.2 mm peaks detected 
by Hill et al. (2005) with SIMBA on the SEST, the factor of 2 gain in angular resolution 
provided by P-ArTeMiS and the absence of confusion by free-free emission (see below) 
allowed us to resolve at least two of the SIMBA peaks (SI and S3) into multiple 450 /im 
dust continuum sources (see online Fig. [2]). 

To further compare our P-ArTeMiS 450 /J,m mosaic with the SIMBA 1.2 mm map 
obtained by HiU et al. (2005), we smoothed the P-ArTeMiS data to the 22" resolution 
of the SIMBA data, regridded the maps onto the same grid, and produced a map of the 
450 /im to 1.2 mm flux density ratio over the region where significant 1.2 mm continuum 
emission was detected by SIMBA. The resulting ratio map is shown in online Fig. [31 Over 
most of the area with significant emission at both 450 /im and 1.2 mm, the S'lfo^m/'^i 2inm 
ratio ranges between ^ 20 and ~ 40. These values are consistent with the Si, oc B^{Td) 
spectral dependence expected for optically thin emission from dust with an emissivity index 
f3 — 2 (e.g. Hildebrand 1983) and a temperature Td ~ 10-30 K. In the eastern part of the 
S3 clump, however, the observed <54|Q^j^/S'p2min ratio is much lower, with values as small 
as ~ 1.5-7.5 (cf. online Fig. [3]). Such values of the flux ratio are inconsistent with dust 
emission unless the dust temperature is unrealistically low, Td ^ 5-7 K, for a parsec-scale 
clump. This reflects the fact that the 1.2 mm emission observed with SIMBA toward the S3 
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Fig. 1. Total-power 450 fim dust continuum mosaic of the massive star-forming region 
NGC 3576 taken with P-ArTeMiS on APEX. The rms in this mosaic is - 0.45 Jy/10"- 
beam. The peak is ^ 47 Jy/10"-beam. Contour levels are 2, 4.5, 8, 13, 25, and 37.5 Jy/10"- 
beam. The color scale shown on the right is in Jy/10"-beam. The grid shown on the left 
corresponds to the instantaneous field of view of the 16 x 16 -pixel P-ArTeMiS array. The 
10" (FWHM) beam size is also indicated. 

clump is strongly contaminated by free-free radiation from the RCW 57 HII region (marked 
by white contours in Fig. [3]). Extrapolating the 1.3 cm continmmr emission measured by 
ATCA (Purcell et al. 2008) to 1.2 mm with the cx spectrum of optically thin 

free-free radiation, wc estimate that free-free emission contributes more than 50% of the 
1.2 mm flux densities observed with SIMBA in the eastern part of the S3 clump. 

3.1. Detection of embedded protostars and derived properties 

Two of the seven compact P-ArTeMiS sources listed in online Table [U S3-M4 and S3-M5, 
are closely associated with 6.7 GHz methanol masers (see Fig. [J) and bright MSX sources. 
The four P-ArTeMiS sources Sl-Ml, S1-M2, S4-M6, S5-M8 were not detected by MSX 
but are closely associated with 22 GHz water masers (see Fig. [J) . The seventh P- ArTcMiS 
source, labeled S3-C3 here, also remained undetected by MSX but coincides with a warm 
{Tnh3 ~ 30-40 K) ammonia core recently identified by Purcell et al. (2008) with ATCA 
(and labeled C3 by these authors). These associations indicate that the ~ 0.1-0.4 pc dusty 
dense cores detected with P-ArTeMiS contain at least one protostellar object each. 

While the ~ 10" (HPBW) angular resolution of P-ArTcMiS at 450 ^im is hardly suf- 
ficient to resolve individual condensations/envelopes within the dust cores at the distance 
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of NGC 3576, our data can be used to estimate the masses and bolometric luminosities of 
the dominant protostellar sources as follows: 

Assuming optically thin dust continuum emission at 450 /im, the gas masses M of the 

S d} 

sources were derived from the measured flux densities S'450um using M = isgt™ — 

where K450/im is the dust opacity per unit (gas + dust) mass column density at A = 
450 /Lim, and -8450^111 is the Planck function Bi,{Td) for the dust temperature Td- We 
adopted K45o^m = 0.04 cm^g~^, which is consistent with the dust opacity law avocated 
by Hildebrand (1983) and appropriate in regions of moderately high gas densities {nH2 ^ 
10^ cm""^ - cf. Henning et al. 1995). The mean dust temperature of each core was estimated 
from the observed 450 /im to 1.2 mm flux density ratio (cf. Fig. [J]) assuming optically thin 
dust emission with an emissivity index [3 = 2 (Hildebrand 1983). The results, given in 
Col. 7 of Table [U generally agree with the gas kinetic temperatures derived from the NH3 
observations of Purcell et al. (2008). Because of free- free contamination in the eastern part 
of the S3 clump (see above) , our method provides only lower limits to the dust temperatures 
of S3-M4 and S3-M5. For these two cores, we thus relied on the gas temperature estimates 
of Purcell et al. (2008). The adopted dust temperatures and the corresponding core masses 
are listed in Col. 2 and Col. 3 of online Tabled respectively. Column 4 provides an estimate 
of the envelope mass, Mg„'!JfP'^, within a diameter of 0.03 pc of the dominant protostellar 
object in each core (where 0.03 pc ~ 6000 AU is the typical fragmentation lengthscale 
observed in nearby cluster- forming regions - cf. Peretto et al. 2007). Mg^^^P'^ was obtained 
by first calculating the mass corresponding to the peak flux density in the 10" beam (given 
in Col. 4 of Table [1] and corresponding to a ^ 0.12 pc aperture at 2.4 kpc), and then 
applying a scaling factor of 0.03/0.12 — 1/4 to that mass, assuming a p cx density 
distribution so that the measured mass scales as the aperture size. A mean column density 
and a mean volume density for each core are also given in Col. 8 and Col. 9 of Tabled 

The bolometric luminosity L^oi of the dominant protostellar object in each core was 
estimated by fitting the observed SIMBA, P-ArTeMiS, and MSX flux densities (or upper 
limits) with the large grid of protostellar spectral energy distributions (SEDs) computed 
by Robitaille et al. (2006, 2007), using their online SED fitting tool. Unfortunately, no 
Spitzer-MYPS observations of NGC 3576 are available yet, consequently resulting in poorly- 
sampled SEDs and large uncertainties on L^oi ■ Coupling a SED fitting approach with our 
P-ArTeMiS 450 fMn. results nevertheless allowed us to derive a useful range of possible Li,oi 
values for each protostellar core, which is given in Col. 5 of Tabled) 

4. Discussion and conclusions 

To put our present results in the context of similar observations of lower-mass Class 
and Class I protostellar objects and to facilitate comparison. Fig. d] shows the locations 
of the dominant sources of NGC 3576 in the Menv~Lboi evolutionary diagram for proto- 
stars (cf. Andre & Montmerle 1994, Saraceno et al. 1996, Andre et al. 2000, Molinari et 
al. 2008). The sources were placed in Fig. Sunder the assumption that the most lumi- 
nous object in each core (dominating Lhoi in Table ^ coincides with the column den- 
sity peak of that core (dominating M^^^^"^). Evolutionary tracks, computed assuming 
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Fig. 4. Envelope mass versus bolometric luminosity diagram comparing the locations of the 
protostellar sources of NGC 3576 (filled and open squares with error bars) with the positions 
of low- mass Class I (filled circles) and Class objects (open circles) (see Andre et al. 2000), 
as well as intermediate-mass protostars in NGC 2264 (triangles - Peretto et al. 2007, Maury 
et al. in prep.). Model protostellar evolutionary tracks, computed for various final stellar 
masses (indicated above the L^oi axis) assuming the accretion/luminosity history described 
in the text, are superimposed. Evolution proceeds from the upper left to the lower right 
as indicated by two arrows on each track, plotted when 50% and 90% of the final stellar 
mass has been accreted, respectively. The straight lines show two Menv-Lboi relations 
marking the conceptual border zone between the Class {M^nv > M^^jt) and the Class I 
{Menv < Mi,/€) stage (where e is the local star formation efficiency). The dashed line is 
such that Menv oc Lboi (cf. Andre & Montmerle 1994), while the dotted relation follows 
Menv OC L'^fi as suggested by the accretion scenario adopted in the tracks. 

each protostar forms from a bounded condensation of finite initial mass Menv{0) and has 
Lboi — GMi,{t)Macc{i) / R*[t) + Li,[t) , where R^{Mi,) is the protostellar radius and L^, is the 
interior stellar luminosity (Stabler 1988, Hosokawa & Omukai 2008), are shown. The mass 
accretion rate and the envelope mass were assumed to be related by Macc{t) — <^Menvit)/T, 
where e = 50% is the typical star formation efficiency for individual cores (cf. Matzner & 
McKee 2000) and r = 10^ yr is the characteristic timescale of protostellar evolution, lead- 
ing to Macc{t) and Menv{i) functions declining exponentially with time (see Bontcmps et 
al. 1996). It can be seen that the two tracks corresponding to final stellar masses of 15 Mq 
and 50 M0, respectively, bracket the locations of the NGC 3576 objects in the M^nv-Lboi 
diagram, suggesting that these objects will evolve into stars of masses between ^ 15 
and ^ 50 M0. The relative positions of the NGC 3576 sources in Fig. [3] also support the 
view that Sl-Ml, S1-M2, S4-M6, and S5-M8 are observed at significantly earlier evolu- 
tionary stages than are S3-M4 and S3-M5, while S3-C3 appears to be at an intermediate 
evolutionary stage. Despite relatively large uncertainties, the values derived for both the 
Menv/L^oi evolutionary indicator (Andre et al. 2000, Bontcmps et al. 1996) and the T^oi 
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indicator (Myers et al. 1998) (see Col. 6 and Col. 7 of Tabled respectively) indeed suggest 
that the cores Sl-Ml, S1-M2, S4-M6, S5-M8 contain at least one massive Class object 
each and that the cores S3-M4, S3-M5 are dominated by a massive Class I object. Recall 
that bona-fide low-mass Class objects are characterized by Menv/L^'oi ^ 0-4 Mq/Lq^ and 
Tboi < 70 K (cf. Andre et al. 2000). 

Our P-ArTeMiS 450 /zm results in NGC 3576 illustrate the unique potential of the future 
ArTeMiS bolometer camera for studies of star-forming regions in the southern hemisphere 
with APEX. With a factor of 9 more pixels than P- ArTeMiS at 450 /im, and significantly 
better detectors, ArTeMiS will be nearly two orders of magnitude faster than P- ArTeMiS 
for mapping extended fields. Furhermore, in addition to 450 /im images, ArTeMiS will 
simultaneously deliver 350 /im and 200 fim images, which will provide much better sam- 
pling of protostellar SEDs near their peak, hence improved bolometric luminosities. The 
colors measured between 450 /im and 200 /im will also lead to improved dust temperature 
estimates and will greatly help to discriminate between self-luminous accreting protostars 
and externally heated, starless cores. Combining extensive LABOCA surveys at 870 /im 
(such as ATLASGAL - SchuUer, Menten et al., in prep.) with follow-up imaging at 200- 
450 /im with ArTeMiS will make it possible to identify and characterize a large number 
of candidate massive protostars in the southern sky. These candidate protostars can then 
be studied in greater detail with ALMA in order to shed light on the formation process of 
high-mass stars. In this way, APEX can truly fulfill its role of "ALMA pathfinder" . 
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Fig. 2. Comparison with the SIMBA/SEST 1.2mm contimmm map of NGC 3576 (Hill et 
al. 2005), overlaid as contours on the P-ArTeMiS 450 /im color image. The white SIMBA 
contours go from 0.3 to 1.5 by 0.3 Jy/22"-beam, and the black SIMBA contours are 3, 
4.5, 9, and 13.5 Jy/22"-beam. Blue stars mark the positions of two 6.7 GHz methanol 
maser sites (Caswell 2004), while red stars mark 22 GHz water maser sites (Purcell et al. 
2008). Note that at least two of the SIMBA peaks (SI and S3) are resolved into multiple 
P-ArTeMiS sources. 




Right Ascension (J2000) 



Fig. 3. Color-scale image showing the variation of the 450 /im to 1.2 mm flux density ratio 
across the NGC 3576 complex. The range of observed intensity ratios is indicated by the 
color scale on the right. The (FWHM) angular resolution of this image is the same as that 
of the SIMBA map (Hill et al. 2005), i.e., - 22". The P-ArTeMiS 450-/im map shown 
in Fig. [T] is superimposed as black contours. The white contours outline the extent of the 
RCW 57 HII region as imaged at 1.3 cm with ATCA (Purcell et al. 2008). 
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Table 1. Protostellar sources detected with P-ArTeMiS in NGC 3576 
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-61° 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Sl-Ml 


±Z _LU.U 


1 ( O i. 


15. ±0.9 


45 ± 2.5 


34 ±2 


21 ± 3 


H2O maser, no MSX 


S1-M2 


12™14.2= 


17'26" 


18. ±0.7 


53 ±2 


27 ±2 


16 ±2 


H2O maser, no MSX 


S3-M4 


ll'"53.2= 


18'27" 


51. ±0.5 


169 ±1.5 


> 20" 


Si 12 


CH3OH maser, MSX 


S3-M5 


11™49.65" 


18'52" 


50. ±0.5 


199 ± 2 


^ 25" 


^ 15 


CH3OH maser, MSX 


S3-C3 


11™47.5" 


18'49" 


45. ±0.5 


192 ± 2 


40 ±2 


33 ±4 


NH3 peali (C3), no MSX 


S4-M6 


11'^ 40.5= 


19'43" 


15.5 ±0.5 


51 ±1.5 


31 ± 1 


19 ± 1 


H2O maser, no MSX 


S5-M8 


ir"34.3= 


21'17" 


11. ±0.5 


33 ± 1 


23 ±2 


13.5 ± 2 


H2O maser, no MSX 



" Lower limit to for the dust source after correcting 51.2111111 for free-free emission. 



Table 2. Derived properties of the NGC 3576 protostellar sources 



Core 
(1) 


Td 
(K) 
(2) 


j^0.35pc 

(Mo) 
(3) 


j^O.OSpc 

(Mq) 
(4) 


Lboi 
(Le) 
(5) 


(6) 


Tbol 
(K) 
(7) 


< Nh2 >0.35pc 

(102=* cm-2) 
(8) 


< riH2 >0.35pc 

(10^ cm~^) 
(9) 


Sl-Ml 


21 


250 


21 


500-10* 


0.1-0.5 


^ 80 


1.4 


2.0 


S1-M2 


16 


540 


45 


500-10* 


0.2-1.1 


^ 70 


3.0 


4.2 


S3-M4 


35 


400 


30 


0.9-3 xlO^ 


0.01-0.03 


120-170 


2.2 


3.1 


S3-M5 


35 


460 


29 


0.5-2.5 xlO'^ 


0.02-0.04 


110-160 


2.6 


3.6 


S3-C3 


33 


490 


29 


0.5-10 xlO" 


0.03-0.2 


^ 60 


2.7 


3.8 


S4-M6 


19 


350 


26 


500-10* 


0.1-0.6 


^ 80 


2.0 


2.7 


S5-M8 


13.5 


500 


42 


500-10* 


0.15-1.0 


^ 90 


2.9 


4.0 



Notes: Column (3) gives the core mass derived from Slso^ni integrated over a 30" diameter 
(Col. 5 of Table [1} ; uncertainty is a factor of ~ 2 on either side due to the uncertain value of 
K450(jm. Column (4) is the envelope mass inferred for the dominant object in each core (see 
text); estimated uncertainty is a factor of ~ 4 on either side. Column (5) gives the bolometric 
luminosity estimated from Sfsojjnii 'S'i'2mm) ^nd MSX flux densities (or upper limits), using the 
fitting routine of Robitaille et al. (2007) . 



